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Abstract 

We define a new scaling limit of matrix models which can be related to the 
method of causal dynamical triangulations (CDT) used when investigating two- 
dimensional quantum gravity. Surprisingly, the new scaling limit of the matrix 
models is also a matrix model, thus explaining why the recently developed CDT 
continuum string field theory (jarXiv: 0802. 071 9t) has a matrix-model representation 
(laDCiv:0804.0252i) . 



1 



1 Introduction 



The great versatility of matrix models or matrix integrals in theoretical physics is 
well illustrated by their particularly beautiful application in two-dimensional Eu- 
clidean quantum gravity (see [11 [21 [Sj H] for reviews). This theory can be defined 
as a suitable sum over triangulations, so-called "dynamical triangulations" (DT), 
whose continuum limit is obtained by taking the side lengths a of the triangles to 
zero. The method of DT was originally introduced as a nonperturbative world- 
sheet regularization of the Polyakov bosonic string O El [7] . There it was used with 
success (or to disappointment, depending on one's taste) to show rigorously that 
a tachyon-free version of Polyakov's bosonic string theory does not exist in target 
space dimensions d > 1 However, when viewed as a theory of 2d quantum 
gravity coupled to matter with central charge c < 1, the theory - noncritical string 
theory - is perfectly consistent, and matrix models have been used to solve in an 
elegant way the combinatorial aspects of the DT construction, where one sums over 
random surfaces glued together from equilateral triangles. 

The DT approach possesses a well-defined cut-off, the length a of the lattice links. 
As has been discussed in many reviews (for instance the ones mentioned above), a 
continuum limit can be defined when the lattice spacing is taken to zero while simul- 
taneously renormalizing the bare cosmological constant and possibly other matter 
coupling constants. However, the continuum limit in question has some unconven- 
tional properties. In this article we will show that there is another way of taking 
the scaling limit of the matrix models, which still relates them to a summation over 
triangulated random surfaces, the so-called causal dynamical triangulations (CDT) 
[9]. We will show that this limit is in a way more natural and does not lead to 
the somewhat unconventional renormalization encountered in the standard DT ap- 
proach. 

In the next section, we will review briefly the conventional situation in the sim- 
plest case, that of pure Euclidean two-dimensional quantum gravity without matter. 
We follow the notations in [3]. In Sec. [3] we introduce a new limit of the correspond- 
ing matrix model and analyze its interpretation in terms of random surfaces. Sec. Hj 
discusses the results obtained and outlines possible applications. 

2 The old matrix model and its continuum limit 

Our starting point is the Hermitian matrix integral 

^(^) = /#e^-^a^M^^)=gl Jd^Pe-'^^-i^') [ftr^^y, (1) 

where is an x Hermitian matrix. This integral is formal since (assuming 
^ > 0) it is not convergent unless an analytic continuation is performed. However, 
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by power-expanding exp{Ng tr0^) as indicated in eq. ([T]) the matrix integral defines 
a formal power series in the coupling constant g. Furthermore, this expansion can 
be given the geometric interpretation of gluing together triangles in all possible 
ways. The size N of the matrix acts as a factor, organizing the power series further 
into a summation over surfaces of fixed topologyo It turns out that the power 
series corresponding to a given fixed topology is convergent with convergence radius 
cjc independent of the topology considered, a property that has made the matrix 
integrals useful in the study of noncritical strings. The particular rearrangement 
of the power series ([1]) according to topology, i.e. in powers of 1/iV^ is called the 
large- expansion. 

Let us briefly review how a "continuum limit" of the matrix integrals can be 
associated with noncritical string theory in the simplest case of c = 0, "pure" Eu- 
clidean two-dimensional quantum gravity, where there is no extended target space 
(equivalently, no matter fields coupled to 2d gravity). This is best done by calculat- 
ing an "observable", the so-called disc amplitude or Hart le- Hawking wave function 
w{z) of the 2d universe to leading order in the expansion. One has 

, ^ 1 1 \ 1 >^ (tr0") 
w{z) = ( — tr \ - \ \ ^ / 



N z-d)/ N ^ z 

n=0 



1 f> (tr0") 



where 

(tr^'^) = Z{g)-' I d<j) tr^ e-^*'"(^'^'-t<^'), (3) 



again to be viewed as a formal power series in g. To leading order in one finds 

w{z) = ^(^z-gz^ + g{z - h)^{z - c){z - d)^ , c> 0, c > (4) 

where the constants h{g), c{g) and d{g) are functions of g and uniquely determined 
by the requirement that w{z) fall off like 1/z (with coefficient 1) for z — > cxo. 

The geometric interpretation of w{z) is as follows: the term in ([2]) corresponding 
to (tr 0") represents the summation over all triangulations with a boundary consist- 
ing of n links and (to leading order in the large- expansion) the topology of a discH 
Consequently, zw{z) represents the sum over all triangulations with the topology 
of a disc and the constant \nz has the interpretation of a boundary cosmological 
constant lb, such that a boundary of length n is assigned a weight e"'*"". In the 
same way as w{z) has a power expansion in g, it also has a power expansion in 1/z 



^By expanding the matrix integral in a power series in g and subsequently in powers of l/N"^, 
we are of course leaving out any nonperturbative contributions to the matrix integral (defined in 
some way, e.g. by analytic continuation) which are not captured by such an expansion. 

^When performing the Gaussian integral, each boundary link, represented by a factor of cj) in 
tr0", is glued either to a triangle or to another boundary link. In the latter case we call it a 
"double link" and it will not be glued to other triangles, see Fig. [T] 
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Figure 1: Graphical representation of relation ([5]): The boundary contains one marked 
link which is part of a triangle or a double link. Associated to each triangle is a weight 
and to each double link a weight 1. 

and the radius of convergence is c{g). The geometric interpretation can be further 
elaborated upon by noting that w{z) satisfies the combinatorial equation 

w{z) = zgw{z) + -w'^{z) + -Q{z,g), (5) 
z z 

where Q{z^ g) is a polynomial in z which is uniquely determinecifl by the requirement 
that the solution to ([5]) have the form (jlj) and fall off like 1/z. The recursion relation 
leading to ([5]) is represented graphically in Fig. [1] and will be important in the next 
section. 

The continuum limit of the matrix model, relevant for noncritical string theory, 
is obtained as follows. Note first that the radius of convergence gc is determined by 
the condition 

bi~g) = cig), (6) 

namely, the point where w{z) changes analytic structure from (z — cY^'^ to (z — cY^"^. 
Next we fine-tune g to gc and z to c{gc) according to 

g = gc{l-Aa^ + 0{a^)), z = c{gc) + aZ + 0{a^). (7) 

These assignments can be viewed as additive renormalizations of the bare cosmo- 
logical and boundary cosmological constants such that A and Z now represent the 
renormalized coupling constants in the limit where the lattice spacing a ^ 0. The 
rationale behind interpreting A as a continuum cosmological constant comes from 
considering the term in w{z) with k triangles glued together, for very large k. It 
will appear as 

^ s{k)e-^'''^, (8) 

■^For a combinatorial interpretation of Q{z,g) see, for instance, [4]. Q{z) is present in eq. ([5]) 
because Fig. [T] is not correct wtien ttie boundary consists of less than two links. 
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where s{k) is a subleading term and we associate A{k) = ka? with a macroscopic 
area. The couphng constant multiplying this area is by definition proportional to the 
continuum renormalized cosmological constant. Similarly one is led to the conclusion 
that Z can be interpreted as a continuum boundary cosmological constant. 
Inserting ([7j) into (jl]), we obtain in the limit of a ^ that 

w{z) = ]^{z-~gz^ + ~gVcCgc)-dCgc) a'^'WEiZ, A) (l + 0(a))) , (9) 

where 

WEiZ,A) = {Z-^/2A/3)^Z + 2^/2A/3 (10) 

is called the continuum Hartle-Hawking wave function or (using string terminology) 
the continuum disc amplitude. 

The above arguments can be generalized from the matrix potential in ([1]) to a 
general potential of the form 

1 " / 

V{<P) = -<p'-gJ2-^\ (11) 

^ ■ 1 ^ 

1=1 

again leading to a formal power series after expanding the exponential in powers 
of the coupling constant g and performing the remaining Gaussian integrals. The 
geometric interpretation is analogous to that for the cubic potential, only that now 
we glue together i-gons and j-gons, where i and j run from 1 to n, with each z-gon 
assigned a relative weight tj. By allowing such polygons we are clearly considering a 
generalization of simplicial complexes, but we will for convenience continue to refer 
to these as "triangulations" . The Hartle-Hawking wave function w{z) now satisfies 
the combinatorial equation 

w{z) = 9(^2 '^^^^''^ ^(^) + i ^'(^) + i (12) 

where again the polynomial Q{z,g) is determined by the requirement that w{z) fall 
off like 1/z and have a single cut. The equation has the graphical representation 
shown in Fig. [2J 

As long as the weights ti are positive (a natural requirement if one wants to 
assign an area to the i-gon), one obtains 

w{z) = i (V\z) + gtWc{gc)-d{g,) a'/^WE{Z, A)) (13) 

in the limit a 0. In (1131) . V'{z) denotes the derivative with respect to z of the 
potential V{z) defined in (ITTl) . g^ the critical value of g (for fixed t^'s), i.e. the radius 
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Figure 2: Graphical representation of relation (jl2p : the marked link of the boundary 
belongs either to an i-gon (associated weight gti) or a double link (associated weight 1). 
It is also a graphical representation of eq. ([38j) if instead of weight 1 we associate a weight 
gs to the marked double link. 

of convergence of the formal power expansion in g for fixed topology (in this case 
that of the disc), and one has made an ansatz similar to (JTj), namely, 

g = gcil-Aa^), z = c{gc) + aZ (14) 

for the coupling constants. One peculiar aspect of ([9]) and (|T3l) is the non-scaling part 
V'{z)/2. This term clearly dominates when a — >■ 0, and in fact renders the average 
number of z-gons present in the ensemble with partition function w{z) finite, even 
at the critical point f|T4l) . This somewhat embarrassing fact can be circumvented 
by differentiating w{z) a sufficient number of times with respect to g and z, which 
will make these "non-universal" contributions vanish. For example, differentiating 
w{z) twice with respect to g yields an expression which diverges in the limit a — > 0, 
allowing one to ignore the finite contribution from non-universal terms, 

dMz) _l 1 1 d^WEiZ,A) 

dg' aV2 ^(^+^/2a73)3/2 a5/2 dA' ' ^ ' 

In this sense one may view We{Z,A) as the continuum disc function. However, it 
would be incorrect to say that w{z) was the regularized Hartle-Hawking wave func- 
tion from which one obtained the continuum Hartle-Hawking wave function through 
an additive renormalization of the cosmological and the boundary cosmological con- 
stants. 

In the following section, we will discuss a new kind of scaling limit for matrix 
models, where such non- universal term are not present. 
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3 The new scaling limit of matrix models 



Hermitian matrix models are often analyzed in terms of the dynamics of their eigen- 
values. Since the action (fTTj) is invariant under the transformation —* U(f)W , with 
U G U{N) a unitary x A^-matrix, one can integrate out the "angular" degrees of 
freedom. What is left is an integration over the eigenvalues k of only, 

i=l k<l 

where the last factor, the Vandermonde determinant, comes from integrating over 
the angular variables, and where 

N 

iiV{ct>) = Y,V{k)- (17) 

Naively one might expect that the large- limit is dominated by a saddle-point with 
V'{1) = 0. However, this is not the case since the Vandermonde determinant in flTB]) 
contributes in the large- limit. The cut which appears in w{z) is a direct result 
of the presence of the Vandermonde determinant. In this way one can say that the 
dynamics of the eigenvalues is "non-classical", deviating from V'{1) = 0, the size 
of the cut being a measure of this non-classicality. We will now introduce a new 
coupling constant Qs in the matrix model by substituting 

Vi<P) ^ - V{<P), (18) 
9s 

and consider the limit Qs 0. As will become clear in due course, this controls 
and reduces the size of the cut and thus brings the system closer to a "classical" 
behaviour. 

In the analysis it will be convenient to keep the coupling constant ti > 0, which 
can be motivated as follows. Consider a "triangulation" consisting of Ti one-gons, T2 
two-gons, T3 triangles, T4 squares etc., up to T„ n-gons, and expand the (^-dependent 
part of the associated matrix integral 



JV 



Z{g,gs) = J d(t)e'Ts''{y'-a^^-^^^l (19) 

in powers of g. Each i-gon appearing in the triangulation has a factor g/ gs associated 
with it. At the same time, the Gaussian integration will produce a factor g^ , where 

L = 1ti + T2 + + ■ ■ ■ + (20) 
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is the number of links in the triangulation, since each Gaussian integration corre- 
sponds precisely to a gluing of an i-gon and a j-gon. The total coupling-constant 
factor associated with the triangulation is therefore given by 

^Ti+--+T„^-Ti/2+T3/2+--+{n/2-l)T„ j^21) 

We observe that in the limit (7^ ^ 0, a necessary condition for obtaining a finite 
critical value gdgs) is Ti > 0. We should emphasize that the analysis described 
below can be carried out also if we suppress the appearance of any one-gons (by 
setting ti = 0) in our triangulations, but it is slightly more cumbersome since then 
9c{gs) — > C)0 as — >• 0, requiring further rescalings. 

For simplicity we will consider the simplest nontrivial model with potential 



V{<l>) = - ^-g<P + -0^ - J (23) 
and analyze its behaviour in the limit gs ~^ 0- The disc amplitude now has the form 

w{z) = ^ [-g + z ~ gz" + g{z - h)^J{z - c){z - d)) , (24) 

and the constants 6, c and d are determined by the requirement that w{z) 1/ z for 
z oc. Compared with the analysis of the previous section, the algebraic condition 
fixing the coefficient of to be unity will now enforce a completely different scaling 
behaviour as (?s — 0. 

For the time being, we will think of gs as small and fixed, and perform the scaling 
analysis for gdgs)- As already mentioned in eq. (EI) above, the critical point gc is 
determined by the additional requirement that b{gc) = c{gc), which presently leads 
to the equation 

{1 - 4g'f^' = I2V3 g', g,. (25) 

Anticipating that we will be interested in the limit gs ^ 0, we write the critical 
points as 

gM = ^(1 - Ag^gs)), AgM = lo'J' + 0{g'J'), (26) 

and 



zM = c{g., gs) = ^ + 1=1+ + Oig^/% (27) 



^In accordance with the remark above this model can be related to the model without one-gons 
by the field and coupling-constant redefinitions 

^=yT^L+- ^=—], ~9 = ^J—. (22) 
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while the size of the cut in ( 123|) . c{gc) — d{gc), behaves as 

cig,)-dig,) = Agl/' + Oig'J'). (28) 

Thus the cut shrinks to zero as — > 0. 

Expanding around the critical point given by fl26l) - fl27j) a nontrivial limit can be 
obtained if we insist that in the limit a ^ 0, gs scales according to 

gs = Gsa\ (29) 

where a is the lattice cut-off introduced earlier. With this scaling the size of the 
cut scales to zero as 4aGy^. In addition {z — c){z — d) oc a if we introduce the 
standard identification z = c{gc) + aZ. This scaling is different from the 

conventional scaling in Euclidean quantum gravity where ^y(z — c){z — d) oc a^/^ 
since in that case {z — c) scales while {z — d) does not scale. 
We can now write 

g = g,{g,){l - a^A) = ^(1 - a^A.^t + 0(a^)), (30) 

with the identifications 

Acdt^A + ^Gf, -g = \, (31) 

as well as 

z = Zc + aZ = z + aZcdt + O(a^), (32) 

with the identifications 

Zcdt = Z + Gl'\ z = l. (33) 
Using these definitions one computes in the limit a — that 



In f lM|) . the constant H (or rather, its rescaled version h = H/^/2K^t ) satisfies the 
third-order equation 

^'-'"^P^Fi-O. (35) 

which follows from the consistency equations for the constants 6, c and d in the limit 
a — * 0. We thus define 

w{z) = -iyedt(^cdt,Acdt,a) = -W{Z,A,Gs) (36) 
a a 

in terms of the continuum Hartle-Hawking wave functions WcAt{ZcAt) ^cAt^Gg) and 
W{Z,KG,). 
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Notice that while the cut of {z — c){z — d) goes to zero as the lattice spacing a, 
it nevertheless survives in the scaling limit when expressed in terms of renormalized 
"continuum" variables, as is clear from eq. Only in the limit — it 

disappears and we have 

«^(^) = -l^cdt(^cdt,Acdt,a) -^^ -— — (37) 
a Gs^o a Zcdt + v2Acdt 

which is the original CDT disk amplitude introduced in [9]. 

A number of comments are in order to put this result into the context of previous 
work concerning the standard and generalized 2d CDT theory, as well as Euclidean 
quantum gravity: 

(1) Eqs. (134|) and (135!) are precisely the equations that were derived in the gener- 
alized CDT model of references [101 EH 112] • This suggests that Gg should be 
interpreted as a coupling constant associated with the splitting of a spatial uni- 
verse into two. This is a meaningful statement in a universe with Lorentzian 
signature, where (in the simplest case) such a splitting is associated with an 
isolated point where the metric and its associated light-cone structure are 
degenerate, which has a diffeomorphism-invariant meaning. This was the mo- 
tivation for considering such processes in CDT in the first place, leading to 
a generalization of the original CDT disc amplitude alluded to in ( 1371) to the 
expression Wcdt(-^cdt, Acdt, G^) given by eq. (IMl) when Gs > 0. However, we 
can also formally make such an association in this purely Euclidean matrix 
model by noting that by introducing the coupling constant Qs (and assuming 
i^i = ^3 = 1 and ti = otherwise) eq. (IT^ is changed to 



wiz) 



] + + lQiz,9)- (38) 



Going back to Fig. 12, this suggests that one should associate a factor instead 
of a factor 1 with the graph with the double line. Geometrically this can 
be viewed as a process where a triangle is removed at a marked link (and 
a new link is marked at the new boundary), except in the case where the 
marked link does not belong to a triangle, but is part of a double-link, in 
which case the double link is removed and the triangulation is separated into 
two. If one thinks of the recursion process in Fig. 12] as a "peeling away" 
of the triangulation as proper time advances, the presence of a double link 
represents the "acausal" splitting point beyond which the triangulation splits 
into two discs with two separate boundary components (i.e. two separate 
one-dimensional spatial universes). Associating an explicit weight Qs with this 
situation and letting gs suppresses this process compared to processes 
where we simply remove an i-gon from the triangulation. Nevertheless, in the 
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limit a — the process survives precisely when we scale (7^ — as prescribed 
by eq. (12^ . The interpretation of this process, advocated in [13], is that it 
represents a split of the spatial boundary with respect to (Euchdean) proper 
time. 

(2) Using eq. (!35|) we can expand w{z) into a power series in Gg/ Aj^^^ whose radius 
of convergence is 1/3 a/3. For fixed values of Acdt, this value corresponds to the 
largest value of Gs where ( l35ll has a positive solution for h. The existence of 
such a bound on Gg for fixed A^dt was already observed in |10j . This bound can 
be re-expressed more transparently in the present Euclidean context, where it 
is more natural to keep the "Euclidean" cosmological constant A fixed, rather 
than Acdt- We have 

Gg 1 3Gs^ / ^ 

(2Acdt)^/2 - ^ ^ 2A + 3Gr - ^' ^^^^ 

which for fixed A > is obviously satisfied for all positive Gg- In order to 
see that the usual Euclidean 2d quantum gravity (characterized by some finite 
value for gg) can be rederived from the disc amplitude flMl) . let us expand eq. 
for large Gg. The square root part becomes 

a-' Gf'^ - v/2A73) ^Jz + 2^2KJ^, (40) 

which coincides with the generic expression a^/'^WE{Z^ K) in Euclidean 2d 
quantum gravity (c.f. eqs. ([9]) and (ITOl) ) if we take Gg to infinity as Qg/a?. 
However, if we reintroduce the same scaling in the \^'(z)-part of 'w{z), it does 
not scale with a but simply goes to a constant. This term would dominate 
w{z) in the hmit a — > if one did not remove it by hand, as is usually done 
in the Euclidean model (see Sec. [2]). 

(3) Why does the potential V'{z) (and therefore the entire disc amplitude w{z)) 
scale (like 1/a) in the new continuum limit with Qg = Gga^, a — 0, contrary 
to the situation in ordinary Euclidean quantum gravity? This is most clearly 
seen by looking again at the definitions fl30|) and fl32l) . Because of the vanishing 

V'iz,g) = 0, V"iz,g) = (41) 

in the point {z,g) = (1, 1/2), expanding around {z,g) according to (l30l) . (!32|) 
leads automatically to a potential which is of order when expressed in terms 
of the renormalized constants (Zcdt,Acdt), precisely like the square- root term 
when expressed in terms of (Zcdt, Acdt)- 

The point {z, g) differs from the critical point {zc{gg), gdgs)), as long as gg 7^ 0. 
In fact, both 1/z and g lie beyond the radii of convergence of l/z and g, which 
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are precisely l/zc^gs) and gdgs)- However, since the differences are of order a 
and a^, respectively, they simply amount to shifts in the renormalized variables, 
as made explicit in eqs. (PT]) and ([33D- Therefore, re-expressing W{Z,A,Gs) 
in fl36|) in terms of the variables Zcdt and Acdt simply leads to the expression 
W^cdt(^cdt, Acdt, Cs), first derived in [10]. Similarly, any geometric quantities 
defined with respect to Z and A can equally well be expressed in terms of Zcdt 
and Acdt- For instance, the average continuum length of the boundary and the 
average continuum area of a triangulation are given by 

- dz - dz:;, ' (^^^ 

_ dlnW{Z,A,Gs) _ 91nWcdt(^cdt,A,dt,a) 
" dA - dl^t • ^^^^ 

In the limit of Gs —>■ 0, the variables {Z,A) and (Zcdt, Acdt) become identical 
and the disc amplitude becomes the original CDT amplitude (Zcdt + V2Acdt)~"^ 
alluded to in ([3] 



(4) In the matrix potential fl2^ . which formed the starting point of our new scaling 
analysis, we are still free to perform a change of variables. Inspired by relations 
(I30|) - (l33|) . let us transform to new "CDT" -variables 



0^^/ + a$ + O(a2), (44) 

at the same time re-expressing g as 

g = g{l-a^Acdt + Oia^)), (45) 

following eq. (l30l) . Substituting the variable change into the matrix potential, 
and discarding a 0-independent constant term, one obtains 

v{<p) = y($) = -^^^^^ (46) 

in the limit a — > 0, from which it follows that 

Z{g,g.,) = a""" Z {Acdt, Gs), ^(Acdt, G,) = J d<!> e'^^^^^*). (47) 

The disc amplitude for the potential V{^) is precisely (Zcdt, Acdt, C^), and 
since by definition 

~ a Zcdt-$' ^^^^ 

the first equal sign in eq. ( l36i) follows straightforwardly from the simple alge- 
braic equation (H6|) . We conclude that the continuum generalized CDT-theory 
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is described by the matrix model with potential V{^). While in the present 
article we have analyzed this only to leading order in A^, the proof that the 
generalized CDT-theory is reproduced by the matrix model with potential 
(H6|) to all orders in was already given in [12]. The important new point 
made here is that there really is a regularized lattice theory with a geometric 
interpretation underlying the formal manipulations of eqs. P6l) - fH7j) . where a 
appears merely as a parameter without an obvious geometric interpretation as 
a cut-off. 

(5) We can generalize the discussion above if we drop the restriction that all ti 
should be positive. Then there exist choices such that 

V'{z,g,U) = --- = V^^\z,g,U) = 0, (49) 



where the derivatives are with respect to z. For a potential satisfying (H9l) we 
obtain a nontrivial limit if we introduce the scaling 

gs = Gsa^, g = g{l-a^A:i^), z = z{l - aZ,at). (50) 

In the limit a — » one obtains 

Acdt) = ^"(^-'^O + ^n(^cdO Vj^^cdt - C){Z,.. - D) ^ ^^^^ 

K(^cdt) = Kil - ^^c\, (52) 

where P„(Zcdt) is a polynomial of order n — 1. The coefficients of P„(Zcdt) 
as well as C and D satisfy algebraic equations which are derived from the 
requirement that VFcdt(^cdt) fall off like l/^cdt- In the limit Gs ^ one finds 

W^edt(^cdt, Acdt, Gs) ^7=^, (53) 

^cdt + V ^Acdt 

which coincides with the original CDT result, valid when no splitting of the 
one-dimensional spatial universe is allowed as a function of proper time. 

The scaling of g in ( l50i) is identical to the scaling for the conventional multi- 
critical one-matrix model p3l [T5] (for reviews see [21 E]). Precisely as in that 
case one has a. {n — l)-dimensional set of deformations around the solution 

(EH). 



4 Discussion 

The CDT model has been advocated as a potential candidate for quantum gravity 
in four dimensions [El |17] , the triangulations in this case constructed by gluing to- 
gether four-simplices in such a way that the geometry is causal in the sense defined 
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in GR. If one does not impose any causal constraint these regularized models have 
so far not led to a continuum theory in four dimensions [TB]. It is of interest to 
study to what extent one can lift the causality constraint in the CDT model and 
study the relation to the pure Euclidean model. This was one of the main motiva- 
tions for the present investigation in two dimensions. The two-dimensional models 
are of course only toy models of quantum gravity, but they have the advantage 
that explicit analytic calculations can be performed. In this spirit it was shown in 
[T2] that the loop equations of the matrix model with the cubic potential (H6|) are 
identical to the continuum string-field equations of a generalized CDT model. The 
generalization consisted in the inclusion of topology changes of the spatial universe, 
which because of their causality- violating nature were excluded from the original 2d 
CDT model. However, it left open the issue of how to understand this generalized 
continuum model as the scaling limit of an ordinary matrix model, with the conven- 
tional geometric interpretation of implementing a gluing of polygons of side length 
a, a being the ultraviolet lattice cut-off, and in addition to understand the precise 
relation of this model to "ordinary" two-dimensional Euclidean quantum gravity. In 
this article we have provided such an interpretation and worked out in detail the 
relation to Euclidean quantum gravity. 

This enables us to view the matrix models with a potential of the type (fTTl) as 
representing regularized path integrals over random surfaces, as one would like to 
do in a theory of two-dimensional quantum gravity. What is rather remarkable is 
that starting with such a matrix model and taking the continuum limit a — > 0, while 
scaling (7^ — > as Gga^, one again ends up with a matrix model, namely, the one 
given by (H6i) . now describing the continuum CDT string field theory formulated in 

In the CDT framework, the new coupling constant Gs describes the splitting 
of a spatial hypersurface into two. Such a concept is not topological in a two- 
dimensional sense, but has a diffeomorphism-invariant meaning in spacetimes of 
Lorentzian signature. Well-behaved causal properties of spacetime were the starting 
point of the original CDT path-integral formulation, motivated by arguments going 
back to [19]. Each individual, causal configuration in the regularized CDT path 
integral allowed a rotation to a unique Euclidean spacetime, leading to a modified 
Euclidean path integral with many fewer geometries than the ordinary Euclidean 
gravitational path integral. What we have shown here is that also this restricted 
class of geometries can be captured by starting out with an ordinary matrix model 
(i.e. one intended to describe Euclidean 2d quantum gravity), and then imposing a 
"penalty" in the form of a coupling constant Qs = Gsd^ for the process where the 
boundary splits into two. 

A potentially powerful application of the model described here concerns matter- 
coupled quantum gravity. It has been shown that in the continuum limit the trian- 
gulations appearing in CDT are much more regular than generic triangulations of 
the disc: a generic triangulation of the disc has fractal dimension four [13], while 



14 



those appearing in CDT have (fractal) dimension two [TOl [20]. As a result, when 
one couples the Ising model to two-dimensional Euclidean quantum gravity the crit- 
ical matter exponents will change from the Onsager values to the so-called KPZ 
values. This change can be traced to the fractal nature of the generic triangulations 
on which the Ising spins are placed, using DT as the framework for (regularized) 
two-dimensional Euclidean quantum gravity. By contrast, if one couples the Ising 
model to Lorentzian quantum gravity in the form of CDT triangulations, one ob- 
tains the Onsager values for the critical exponents as has been shown numerically 
[2T] and using a high-T expansion [22] . By generalizing the framework of this article 
from a one-matrix to a two-matrix model description, one may be able to obtain a 
two-matrix model description of CDT coupled to the Ising model and in this way 
prove analytically that the critical exponents of the Ising model on CDT are indeed 
the Onsager exponents. If this can be done, it might be by far the simplest way to 
calculate explicitly the Onsager exponents of an Ising lattice model. Work in this 
direction will be reported elsewhere. 

It is also somewhat surprising that the matrix integral fH7j) . with the same dimen- 
sions associated to the fields and coupling constants, is encountered in the so-called 
Dijkgraaf-Vafa correspondence [23j. From the gauge theory side V{^) is then the 
tree-level superpotential of the adjoint chiral field $, which breaks the supersym- 
metry of the unitary gauge theory from A/" = 2 to A/" = 1. If one demands that 
this tree-level potential corresponds to a renormalizable theory then the form of the 
potential V{^) is essential unique (and given by (l46l)). Gs is a dimension three 
coupling constant coming either from three of the compactified dimensions in topo- 
logical string theory or alternatively, via the DV-correspondence, from the glueball 
superfield condensate in the gauge theory. In this formalism a continuum interpre- 
tation is associated with the equation defining disk amplitude, rather than with the 
(discretized) surfaces defined by the matrix integral, the equation defining an alge- 
braic surface. This point of view is also adopted in [2S1 126j where the calculation 
of matrix integrals like flTTI) has been carried to a new level of perfection. In these ap- 
plications one encounters eventually much more complicated situations with higher 
order potentials V{^) which admit multi-cut solutions. It would be interesting if 
these multi-cut solutions could also find a random surface interpretation. 

Acknowledgment 

JA, RL, WW and SZ acknowledge the support by ENRAGE (European Network 
on Random Geometry), a Marie Curie Research Training Network in the Euro- 
pean Community's Sixth Framework Programme, network contract MRTN-CT- 
2004-005616. RL acknowledges support by the Netherlands Organisation for Sci- 
entific Research (NWO) under their VICI program. 



15 



References 



[1] F. David, Lectures given at Spring School on String Theory and Quantum 
Gravity, Trieste, Italy, Apr 23 - May 1, 1990. Published in Trieste Spring School 
1990:55-90. 



P. H. Ginsparg and G. W. Moore, arXiv:hep-th/9304011 



P. Di Francesco, P. H. Ginsparg and J. Zinn- Justin, Phys. Rept. 254 (1995) 1 
(^rXiv:hep-th/9306153] . 

J. Ambj0rn, B. Durhuus and T. Jonsson, Cambridge, UK: Univ. Pr., 1997, 
( Cambridge Monographs in Mathematical Physics ), 363 p 

J. Ambj0rn, B. Durhuus, J. Frdhhch, Nucl. Phys. B 257 (1985) 433-449; 
J. Ambj0rn, B. Durhuus, J. Frohlich and P. Orland, Nucl. Phys. B 270 (1986) 
457-482. 

F. David, Nucl. Phys. B 257 (1985) 543; A. Billoire and F. David, Nucl. Phys. 
B 275 (1986) 617. 

V. A. Kazakov, A.A. Migdal and I. K. Rostov, Phys. Lett. B 157 (1985) 295. 

J. Ambj0rn and B. Durhuus, Phys. Lett. B 188 (1987) 253-257. 

J. Ambj0rn and R. Loll, Nucl. Phys. B 536 (1998) 407-434 [hep-th/9 805108]. 

J. Ambj0rn, R. Loll, W. Westra and S. Zohren, JHEP 0712 (2007) 017 
[irXiv:0709.2784l [gr-qc]]. 

J. Ambj0rn, R. Loll, Y. Watabiki, W. Westra and S. Zohren, JHEP 0805 (2008) 
032 |arXiv:0802.07T9l [hep-th]]. 

J. Ambj0rn, R. Loll, Y. Watabiki, W. Westra and S. Zohren, Phys. Lett. B 665 
(2008) 252 [arXiv:0804.0252 [hep-th]]. 

H. Kawai, N. Kawamoto, T. Mogami and Y. Watabiki, Phys. Lett. B 306 
(1993) 19-26 |hep-th/9 302133] ; Y. Watabiki, Nucl. Phys. B 441 (1995) 119 
farXiv:hep-th /9401096]; J. Ambj0 rn and Y. Watabiki, Nucl. Phys. B 445 
(1995) 129 |a rXiv:hep-th /9501049|; S. Catterall, G. T horleifss on, M. J. Bow- 
ick and V. John, PhysV^Lett/ B 354 (1995) 58 |arXiv:hep-lat/9504009]; 
J. Ambj0rn, J. Jurkiewicz and Y. Watabiki, Nucl. Phys. B 454 (1995) 313 
|arXiv:hep-lat/9507014|; H. Aoki, H. Kawai, J. Nishimura and A. Tsuchiya, 
Nucl. Phys. B 474 (1996) 512-528 |h ep-th/9511117] . 

[14] V. A. Kazakov, Mod. Phys. Lett. A 4 (1989) 2125. 



16 



[15] M. Staudacher, Nucl. Phys. B 336 (1990) 349. 



[16] J. Ambj0rn, A. Gorlich, J. Jurkiewicz and R. Loll, Phys. Rev. Lett. 100 (2008) 
091304 (arXiv:0712.2 485 [hep-th]]: [i?Xiv:0807.4 481 [hep-th]. 

[17] J. Am bj0m, J. Jurkiewicz and R. Loll, Phys. Rev . D 7 2 (2005) 064014 
(hiph/0505 154|; Phys. Rev. Le tt. 95 (20 05) 171301 [hep-th/ 0505113|; Phys. 
Lett. B 607 (2 005) 205-213 |hep-th/0411152| ; Phys. Rev. Lett. 93 (2004) 131301 
|hep-th/0404 156]. 

[18] P. Bialas, Z. Burda, A. Krzywicki and B. Petersson, Nucl. Phys. B 472 (1996) 
293-308 (hep-lat /9601024j. 

P. Bialas, Z. Burda, B. Petersson and J. Tabaczek, Nucl. Phys. B 495 (1997) 
463-476 [hep-lat/9608030|. 

S. Catte rall, J.B. Kogut and R. Renken, Phys. Lett. B 328 (1994) 277-283 
|hep-lat/94 01026|; Phys. Lett. B 416 (1998) 274-280 [ hep-lat/970 9007]. 
B.V. de Ba kker and J. Smit, Nucl. Phys. B 439 (1995) 239-258 
|hep-lat/9407014l . 

[19] C. Teitelboim, Phys. Rev. Lett. 50 (1983) 705-708; Phys. Rev. D 28 (1983) 
297-309. 

[20] J. Ambj0rn, J. Correia, C. Kristjansen and R. Loll, Phys. Lett. B 475 (2000) 
24 |arXiv:hep-th/9912267j; J. Ambj0rn, R. Janik, W. Westra and S. Zohren, 
PhysTXett. B 641 (2006) 94 |arXiv:gr-qc/0607013| . 

[21] J. Ambj0rn, K.N. Anagnostopoulos and R. Loll, Phys. Rev. D 60 
(1 999) 10403 5 |arXiv :hep-th/9904012 |; Phys. Rev. D 61 (2000) 044010 
|arXiv:hep-lat/9909129|; J. Ambj0rn, K.N. Anagnostopoulos, R. Loll and 
I. Pushkina, arXiv:0806.350 6' [hep-lat]. 

[22] D. Benedetti and R. Loll, Phys ica A 377 (2007) 373 [arXiv:hep-lat/06 03013|; 
Gen. Rel. Grav. 39 (2007) 863 [arXiv:gr-qc/0611075|. ' ^ 



[23] R. Dijkgraaf and C. Vafa, Nucl. Phys. B 644 (2002) 3-20 
(hep-th70 206255|; |h ep-th/0208048l . 



[24] B. Eynard, JHEP 0411 (2004) 031 |arXiv:hep-th/0407261] . 

[25] L. Chekhov and B. Eynard, JHEP 0603 (2006) 014 [arXiv:hep-th/0504116| . 



[26] B. Eynard and N. Orantin, [arXiv:math-ph/0702045| JHEP 0806 (2008) 037 
|arXiv:07 10.0223' [hep-th]]. 



17 



